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Abstract 

The Next to Minimal Supersymmetric Standard Model (NMSSM) with a Higgs- 
singlet coupling A close to unity, moderate tan/3 and stop masses below 1 TeV 
<t^ minimizes the fine tuning of the electroweak scale, making possible that the 

CN) lightest new particles, except perhaps for the LSP, be part of the extended Higgs 

system. We show how the measurements of the couplings of the 126 GeV Higgs 
boson constrain the region of the physical parameters of a generic NMSSM most 
relevant to this context. In the same region of parameter space we determine the 
cross section for the production of a heavier CP-even scalar together with its 
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total width and its most relevant branching ratios. For comparison we show the 
same analysis for the MSSM. We also comment on how a coupling A > 1 can be 
compatible with gauge coupling unification. 
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1 Introduction 

The lack of signals so far from the direct production of supersymmetric particles calls for 
a reconsideration of the strategy to search for supersymmetry. Focussing on naturalness, 
which to us still looks as the best motivated guideline, special attention is paid to the s- 
particles that have the largest influence on the quadratic terms in the Higgs potential, namely 
the two stops, the left-handed bottom, the higgsinos and, indirectly through the stops, the 
gluino [l]44]. 

The mounting evidence for a Higgs boson /ilhc at 126 GeV with Standard-Model-like 
properties |5]49] has also a significant impact on these issues. While the identification of the 
new resonance with the lightest Higgs boson of a supersymmetric model is a well motivated 
possibility, the measured mass is in serious tension with maintaining naturalness in the 
Minimal Supersymmetric Standard Model (MSSM), requiring either large stop masses, above 
a TeV, in the case of negligible mixing, or a large trilinear A t term. At the heart of this 
problem is the fact that the quartic terms of the Higgs potential in the MSSM are controlled 
by the weak gauge couplings. 

The situation is different in the Next to Minimal Supersymmetric Standard Model (NMSSM), 
where a singlet superfield S couples to the Higgs superfields, H u and H d , via the Yukawa-like 



coupling XSH u H d [10] (see 11 for a review). On one hand, the inclusion of this coupling 
allows to get a 126 GeV Higgs boson mass with both the stop masses well below a TeV. On 
the other hand, values of A > 1 suppress the sensitivity of the Higgs VEV with respect to 
changes in the soft supersymmetry breaking masses, thus still keeping the fine tuning at a 
moderate level even for stop masses up to 1 TeV 12-14 HWe shall comment below on how 



a coupling A > 1 can be made compatible with gauge coupling unification. 

In general terms, to see whether the newly found resonance at 126 GeV is part of an 
extended Higgs system is a primary task of the current and future experimental studies. 
Given the above motivations, this appears to be especially true for the extra Higgs states of 
the NMSSM, which might be the lightest new particles of a suitable supersymmetric model, 
except perhaps for the LSP. A particularly important question is how the measurements of 
the couplings of /ilhc, current and foreseen, bear on this issue, especially in comparison with 
the potential of the direct searches of new Higgs statesjj 

Not the least difficulty that one encounters in attacking these problems is the number of 
parameters that enter the Higgs system of the NMSSM, especially if one does not want to 
stick to a particular version of it but rather wishes to consider the general case. Here we aim 
at an analytic understanding of the properties of the Higgs system of the general NMSSM, 
trying to keep under control as much as possible the complications due to the proliferation 
of model parameters and avoid the use of benchmark points. 

The content of the paper is the following. In Section 2 we establish some relations between 
the physical parameters of the CP-even Higgs system valid in the general NMSSM. In Section 
3 and 4 we consider two limiting cases in which one of the CP-even scalars is decoupled, 
determining in each situation the sensitivity of the measurements of the couplings of /ilhc, 
current and foreseen, as well as the production cross sections and the branching ratios for 



X A recent analysis 15 finds that the fine tuning in the NMSSM with a scale invariant superpotential can 
be above 5% for stop masses up to 1.2 TeV and gluino masses up to 3 TeV for A « 1, moderate tan/3 
and the messenger scale at 20 TeV. 
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For recent studies see e.g. 16 17 



the new intermediate scalar. In Section 5 we compare one of these NMSSM cases with the 
much studied MSSM, using as much as possible the same language. In Section 6 we illustrate 
a possible simple and generic extension of the NMSSM that can make it compatible with 
standard gauge unification even for a coupling A > 1. Section 7 contains a summary and our 
conclusions. 

2 Physical parameters of the CP-even Higgs system in 
the general NMSSM 

Assuming a negligibly small violation of CP in the Higgs sector, we take as a starting point 
the form of the squared mass matrix of the neutral CP-even Higgs system in the general 

NMSSM: 



m| cos 2 /3 + m\ sin 2 /3 (2v 2 X 2 — m\ — ra|) cos j3 sin /3 vM\ 

M 2 = \ (2v 2 \ 2 - m 2 A -m 2 z ) cos f3 sin (3 m\ cos 2 /3 + m§ sin 2 (3 + S 2 vM 2 \ (2.1) 



vM 1 vM 2 Ml 



in the basis % = {H®, H®, S) T . In this equation 

m 2 A = m 2 H± -m 2 w + AV, (2.2) 

where m#± is the physical mass of the single charged Higgs boson, v ~ 174 GeV, and 

5 2 = A 2 K 1 ,m t - 2 ,^-)/sin 2 /3 (2.3) 

is the well known effect of the top-stop loop corrections to the quartic coupling of H u , 
with mi and 9i physical stop masses and mixing. We neglect the analogous correction to 



(2.2) 18 , which lowers m#± by less than 3 GeV for stop masses below 1 TeV. We leave 



unspecified the other parameters in (2.1), Mi, M 2 , M 3 , which are not directly related to 
physical masses and depend on the particular NMSSM under consideration. 

The vector of the three physical mass eigenstates % p h is related to the original scalar fields 
by 

H = R^RfRl 3 n ph = RH phj (2.4) 

where R 1 q is the rotation matrix in the ij sector by the angle 9 = a,j, a. 
Defining 7i ph = (h 3 ,hi,h 2 ) T , we have 

i? T M 2 i? = diag(m 2 3 ,m 2 i ,m 2 2 ). (2.5) 

We identify hi with the state found at LHC, so that m^ = 125.7 GeV. For simplicity we shall 
always consider hi as the lightest CP-even state, although other cases may be compatible 
with current dataPJFrom (2.4) hi is related to the original fields by 



hi = c^(-s a H d + c a H u ) + SjS, (2.6) 

where se = sin6, eg = cos8. Similar relations, also involving the angle a, hold for h 2 and /13. 



3 Our analysis can be applied as well to the case of a new state below /ilhc- See 19 and references therein 
for a recent study of the NMSSM that emphasizes especially this case. 
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Figure 1. Fit of the measured signal strengths of h\ 



- ^lhc- Left: 3-parameter fit of tan/?, 5,5 and 
, r ^ — ^~„~~ w b ^. u „,„ „„,^.^. ^ to ..^ ^ „ 7 - (black), 0.15 (dark grey), 0.3 (light grey). 
The regions overlap in part, but their borders are also shown. Right: fit of s 2 in the case of 5 = 0. 



s 2 . The allowed regions at 95%C.L. are given for s 2 



These angles determine the couplings of h\ = /^lhc to the fermions or to vector boson 
pairs, VV = WW, ZZ, normalized to the corresponding couplings of the SM Higgs boson. 
Denning 5 = a — (3 + tc/2, they are given by (see also [20]) 



Qh x tt 

9htt 



c-y{c & + 



s& 



tan/3' 



Qhibb 
9hbb 



Cj(cs — s.jtan/3), 



9hiVV_ 

9hvv 



C/yCjJ. 



(2.7) 



The fit of all ATLAS [7J, CMS [8] and TeVatron [9] data collected so far on the various signal 
strengths of /ilhc gives the bounds on 5 for different fixed values of 7 shown in Figure [T] 
left and the bound on 7 for 5 = shown in Figure [T] right. To make this fit, we adapt the 
code provided by the authors of 21 . As stated below, we do not include in this fit any 



supersymmetric loop effects. Note that in the region of ss close to zero, a larger s 2 forces 5 
to take a larger central value. 



The matrix equation (2.5) restricted to the 12 sector gives three relations between the 
mixing angles a, 7, a and the physical masses rnh 1: h 2 ,h 3 ^ rn H ± f° r anv given value of A, tan/3 



and A t . In terms of the 2x2 submatrix M 2 in the 12 sector of Ai 2 , eq. (2.1 ), these relations 
can be made explicit as 



det M 2 + m 2 hl 1 



m 2 hi -trM 2 ) 



mt 



in 



hi 



)( m hl 



m 



m 



hi 



m 



h 3 J 

det M 2 + m 2 h3 (m 2 h3 



trM 2 



m 



hi 



m 2 h3 detM 2 



m l 2 m l 3 



+ ™l( m L+ m L-^M 2 ) 



(2.8) 
(2.9) 
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sin2« = ( ± 2|s 7 s -|a/1 — s^yl — sin 2 2£ (m 

+ K 3 - m lA + 4(i + 4)( m L - m D - (! - 4) m l] sin2 

x ( [m 2 ^ - m^ + s 2 (m 2 hl - m 2 h2 )] 2 + (m 2 ha - m 2 h2 )(l - s 2 )sl 
x [2m 2 hi (l + s 2 ) - 2(m 2 3 + S 2 m 2 J + S 2 (m 2 3 - m 2 J(l - S 2 )] 



(2.10) 



where £ is the mixing angle of M 2 . In the limit of s 



0, the angle £ goes to a. There 



are two solutions for a given by the ± signs in (2.10). 



These expressions for the mixing angles do not involve the unknown parameters M±, M 2 , M3, 
which depend on the specific NMSSM. Their values in particular cases may limit the range 
of the physical parameters m hlth2jfl3 ,m H ± and a, j,a but cannot affect eqs. (2.8, 2.9, 2.10). 



To our knowledge, analytical expressions for the mixing angles in the general NMSSM have 
not been presented before. 

To simplify the analysis we consider two limiting cases: 



H decoupled: In (2.1 ) m\ 3> vMi, vM 2 or m/j 3 3> mh lt h 2 an d er, 5 = a — (3 + 7r/2 — > 0, 



Singlet decoupled: In (2.1) M3 3> vMi,vM 2 or m/> 2 3> vrih u h 3 an d a, 7 — )■ 0, 



but we use (2.8, 2.9, 2.10) to control the size of the deviations from the limiting cases when 



the heavier mass is lowered. 

When considering the couplings of the CP-even scalars to SM particles, relevant to their 
production and decays, we shall not include any supersymmetric loop effect other than the 



one that gives rise to (2.3 ). This is motivated by the kind of spectrum outlined in the previous 



Section, with all s-particles at their "naturalness limit" , and provides in any event a useful 
well defined reference case. We also do not include any invisible decay of the CP-even scalars, 
e.g. into a pair of neutralinos. To correct for this is straightforward with all branching ratios 
and signal rates that will have to be multiplied by a factor V /(T + T xx ). Finally we do 
not consider in this paper the two neutral CP-odd scalars, since in the general NMSSM 
both their masses and their composition in terms of the original fields depend upon extra 
parameters not related to the masses and the mixings of the CP-even states nor to the mass 
of the charged Higgs. 



3 H decoupled 



To study this limiting case, it is best to go in the basis (H, h, s) with H = —spH^ + cpH u 
and h = cpHd + spHu, and let H m h 3 decouple, so that a, 5 = a — (3 + n/2 — > 0. For the 
remaining non- vanishing angle 7 one has 



sin 2 7 



m lh ~ m l x 



m 



h 2 



mi 



where 



m 



,2 2 



hh 



m ~Z c 2f3 + A f S 2/3 + A 



(3.1) 



(3.2) 



is the first diagonal entry in the square mass matrix of the reduced basis (h, s). 





Figure 2. H decoupled. Isolines of sin 7. Left: A = 0.8. Right: A = 1.4. The colored region is 
excluded at 95%C.L. by the experimental data for the signal strengths of h\ = /ilhc- 



Under the conditions specified in the previous Section it is straightforward to see that 
the couplings of h\ = /ilhc and h% to fermions or to vector boson pairs, VV = WW, ZZ, 
normalized to the same couplings of the SM Higgs boson, are given by 



Qhiff _ ginvv 



« 



9hvv 



U-y, 



9h 2 ff _ 9h 2 VV _ 



« 



9hvv 



= —s 



■ 7 . 



(3.3) 



As a consequence none of the branching ratios of h\ gets modified with respect to the SM 
ones, whereas its production cross sections, or the various signal strengths, are reduced by a 
common factor e£ with respect to the SM ones with m/ lSM = m^ ■ The fit of all experimental 
data collected so far gives the bound on s^ shown in Figure HI right. 

Upon use of (3.1) the impact of this bound on the parameter space is shown in Figure [2 
for A = 0.8 and 1.4, together with the isolines of different values of si that might be probed 
by future improvements in the measurements of the h\ signal strengths. Larger values of A 
already exclude a significant portion of the parameter space at least for moderate tan/3, as 
preferred by naturalness. In this Section we are taking a fixed value of A t = 73 GeV in (2.3), 

45° 
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which is obtained, e.g., for m^ = 600 GeV, m ( - 2 = 750 GeV and mixing angle 9 t 

As long as one stays at A t < 85 GeV, in a range of moderate fine tuning, and A > 0.8, our 

results do not depend significantly on A t . 

In the same (tan/3, m^) plane of Figure^ and for the same values of A, Figure p^ shows the 
gluon-fusion production cross sections of h% at LHC for 8 or 14 TeV c.o.m. energies, where 
we rescaled by (? the NNLL ones provided in 
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. All other hi production cross sections 
relative to the gluon-fusion one, scale as in the SM with rrih SM = rn h2 . 
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Figure 3. H decoupled. Isolines of gluon fusion cross section a(gg — > /12) at LHC8 and LHC14, 
for the values A = 0.8 and A = 1.4. The colored region is excluded at 95%C.L. 





Figure 4. H decoupled. Isolines of BR(/i2 — > hh). Left: A = 0.8 and v$ = 2v. Right: A = 1.4 and 
vs = v. The colored region is excluded at 95%C.L. 





Figure 5. H decoupled. Isolines of BR(/i 2 ->■ W+W~). Left: A = 0.8 and v s = 2v. Right: A = 1.4 
and vs = v. The colored region is excluded at 95%C.L. 



To determine the decay properties of h 2 it is crucial to know its coupling (g h2h 2/2)h 2 hl 
to the lighter state. In the general NMSSM and in the large tjih limit considered in this 
Section, the leading A 2 -term contribution to this coupling, as well as the one to the cubic 
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Figure 6. H decoupled. Isolines of the total width I\ 2 (GeV). Left: A = 0.8 and vs = 2v. Right: 
A = 1.4 and v$ = v. The colored region is excluded at 95%C.L. 



/ii-coupling (g h 3/6)hl, are given by 



A v ( Vs vg 2 \ 

g hoh 2 = — -= 4 — cos 7 + 12 — cos 37 — 7 sin 7 + 12 cos Ap cos 7 sin 7 + 9 sin 37 
1 8\/2 V v v J 



9h\ 



8^ 
3 

V2v 
X 2 v 2 



A 2 cos 2 7 sin 7, 



im 



hi 



( 2 V S \ A 2 

cos 7 10 — 4 cos 4p cos 7 — 6 cos 27 + 8 — sin 27 ) H 5- cos 

V v / mi 



7, 



(3.4) 



(3.5) 



where vs is the VEV of the singlet. Figures |6] and |4] show the total width of hz and its 
branching ratio into a pair of light states for some choices of vs- The other most significant 
decay mode of \\2 is into a P^-pair, with a branching ratio given in Figure [5j Figure [7] shows 
the triple /ii-coupling normalized to the SM one. 

These results depend on the value taken by Vs, in particular we note that the Higgs fit still 
allows the triple Higgs coupling to get a relative enhancement of a factor of a few (with a 
negative or positive sign) with respect to the Standard Model one, thus yielding potentially 



large effects in Higgs pair production cross sections 24 



In the general case, when H is not completely decoupled, the three angles 5, 7 and a can all 
be different from zero, and the three masses m^, m^ and m#± are all virtually independent. 
In Figure M we show the excluded regions in the plane (tan/?, m^) for m hz = 750 GeV and 



A = 1.4, setting s 2 to two different values in order to fix mn±. When s 2 = one recovers the 
previous H decoupled case in the limit rrih 3 —¥ 00. With respect to this case, both 7 and S 
are free parameters in the fit to the couplings of /ilhc, an d as a consequence the bounds are 
milder than what is expected from using only 7. If s 2 7^ 0, h 2 and h 3 are not decoupled, and 
their masses can not be split too much consistently with all the other constraints. This is 
reflected in a broader excluded region for low m^ in Figure fright, where we take s\ = 0.25. 





Figure 7. H decoupled. Isolines of ghhh/9hhh- Left: ^ = 0-8 an d vs = 2v. Right: A = 1.4 and 
vs = v. The colored region is excluded at 95%C.L. 





Figure 8. H "almost decoupled" with A = 1.4 and rrih 3 = 750 GeV. The dashed isolines are for 
m,fj±. Left: sin 2 a = 0. Right: sin 2 a = 0.25. The colored region is excluded at 95%C.L. In the grey 
area there is no solution for 5. The thick line shows the naive exclusion limit from s 2 only. 



10 



; 145 

i / 

V / 


,135 


200 / / 

; / 300 






tan/3 




tan/3 



Figure 9. Singlet decoupled. Left: isolines of rrihh ( |3.2[ ) (solid), the grey region is unphysical due 
to m'fj± < 0. Right: isolines of A in the region with rrihh — 126 GeV. The dashed isolines are for 
m H ±. The colored regions are excluded at 95%C.L. 



4 Singlet decoupled 



Let us consider now the limit in (2.1) M% S> vMi,vM2, which corresponds to vri2 S> m\ r 



and er, 7 — > 0. In this case the three relations that have led to (2.8, 2.9, 2.10) become 



sin 2a = sin 2/3 



2\ 2 v 



2„,2 



m. 



m 



A\m hl 



+ m% + 8? — 2m? ' 



m 



m 



h 3 



m A\m hl +m 7 + 8 



A\m hl 

2 
Z 



In 



m 



iii ' 



where 



m, 



X 2 v 2 (X 2 v 2 — m 2 z ) sin 2f3 — m 2 hi (m hl — m z — 8 2 ) — m 2 z 8 2 cos 2 f3 



'•m h 



m hh ~ m hi 



(4.1) 
(4.2) 

(4.3) 



Note the difference with respect to the single relation (3.1 ) of the previous Section. Identifying 



as before h\ with the resonance found at the LHC, this determines the mass of h 3 (and of 
H ) for any given value of A and tan/3. Still taking A t = 73 GeV, this allows to show similar 
plots to the ones of the previous Section without having to fix the value of A. 
Defining as before 8 = a — f3 + tt/2, the couplings of h 3 become 



9h 3 tt 

9htt 



sin 8 — 



cos 8 
tan/3' 



9h 3 bb 
9hbb 



sin 8 + tan (5 cos 8, 



9h 3 VV 

9hvv 



sin 8. 



(4.4) 



The allowed regions in the plane (tan/3, A) shown in Figure M left are determined by a 2- 
parameter fit of tan/3, sin 8. This fit results in an allowed region which is virtually the same 
as the one with 7 = in Figure [I] left. When inverting A as a function of tan/3, rrih :i , there 
are two solutions. In Figure [9] right, we show only the one which corresponds to the narrow 
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Figure 10. Singlet decoupled. Isolines of gluon fusion production cross section a{gg — > h 3 ). The 
colored regions are excluded at 95%C.L., and the dashed line shows m H ± = 300 GeV. Left: LHC8. 
Right: LHC14. 



allowed region with rrihh close to 126 GeV. Note that A is restricted to relatively small values. 
As a consequence the analysis becomes more sensitive to values of A t at or above 80 GeV. 
The other two allowed regions in Figure M left, when translated to the (tan j3, mh 3 ) plane, 
correspond to the other solution for A, and are not displayed in Figure [9] right. The first 
one with tan/3 ~ 3 and A > 1 always implies a charged Higgs mass m H ± below 150 GeV, 



which is disfavored by indirect constraints 25 . Note that this region, corresponding to the 



allowed region with large 5 in Figure [TJ is mainly allowed because of the large error in the 
measurement of the bb coupling of /ilho Reducing this error down to about 30% around 
9hibb/9hbb — 1 would exclude the region. The second one with tan/3 ~ 1, allowed by the 
Higgs fit, is also excluded by the bounds on m H ±. As already said, we do not consider here 
the region with rrihh < 126 GeV, which requires a CP-even state lighter than /ilhc- 



means of 20 



The couplings (4.4) allow to compute the gluon-fusion production cross section of h 3 by 



a(gg->h 3 ) = a bM (gg^H(m h3 )) 



A t 



9h 3 tt 



A, 



9h 3 



9hbb 



(4.5) 



where 



A 



Fi(rt,b) 



mt 



t.b 



Flirt) 

2 



Ti 



m 



h 3 



and Fi (r) is a one-loop function that can be found e.g. in 

in Figure 10, where we used the values of <r SM at NNLL precision provided in 23 



26 



27 



(4.6) 

This cross section is shown 
and the 



running masses m^ at NLO precision. We checked the validity of this choice by performing 
the same computation both with the use of masses at LO precision and K-factors 28 , and 



with the program HIGLU 29 30 , finding in both cases an excellent agreement. 



9h 3 h 2 9 9h :i s 

The coupling of h 3 to the lighter state — ^-^3^1 and the triple Higgs coupling —r-h l are 
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Figure 11. Singlet decoupled. Left: isolines of the total width r/ l3 (GeV). Right: isolines of BR(/i3 
hh). The colored regions are excluded at 95%C.L., and the dashed line shows m#± = 300 GeV. 
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Figure 12. Singlet decoupled. Left: isolines of BR(/i3 — > tt). Right: isolines of BR(/i3 
colored regions are excluded at 95%C.L., and the dashed line shows m H ± = 300 GeV. 
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Figure 13. MSSM. Left: isolines of rrihh (solid) and m#± (dashed), the gray region is unphysical 
because of m? A < 0. Right: isolines of A 4 (GeV) (solid) and m H ± (dashed). Light colored regions 
are excluded at 95%C.L. by the Higgs fit, the red region is excluded by CMS direct searches for 
A,H ^t + t-. 



Figures 11 and 12 show the most relevant widths of h%. 



5 The MSSM for comparison 

It is instructive to compare the results of the previous Section with the much studied case of 
the MSSM, using the same language as much as possible, since the MSSM is the A = limit 



of the NMSSM in the singlet-decoupled case. A recent paper 31 analyzes the Higgs system 



of the MSSM in a way similar to ours and gives comments about the heavy Higgs searches 
in different channels (see also [32]). 

A first important difference of the MSSM versus the NMSSM is in a minimum value of 
A t > 85 GeV that is needed to accommodate the 126 GeV Higgs boson as the lightest 
CP-even neutral scalar. Also for this reason, and because we have one parameter less than 
in the previous Section, we let A t vary. As a consequence, in analogy with Figure [9j we 



show in Figure 13 the allowed regions by current experimental data on the signal strengths 
of hi = /ilhc- Note that in the plane (tan/3, m^) the isolines of A t are increasingly large 
at lower tan /3: a sign of increasing fine tuning. Finally, in analogy with Figures |lO]fT2"t we 



show in Figures [14 (16 the gluon-fusion production cross sections and the widths of h 3 for the 
MSSM case. For the production cross sections we have adopted the same procedure of the 
Singlet decoupled case, and performed a further check of our results with the ones recently 
presented in 33 and 31 , finding a very good agreement. 

In the MSSM m\ in (4.3) at A = is the squared mass of the neutral pseudoscalar A, 



unlike the case of the general NMSSM, where m A in (4.3) is only an auxiliary quantity. In the 



same (tan/3, m^ 3 ) plane a(gg — >■ A) is therefore also determined, which allows to delimit the 
currently excluded region by the direct searches for A, h 3 — > t + t~ . Such a region is known 
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tan/3 

Figure 14. MSSM. Isolines of gluon fusion production cross section a(gg — > /13). Light colored 
region is excluded at 95%C.L., the red region is excluded by CMS direct searches for A,H — >• t + t~ . 
Left: LHC8. Right: LHC14. 




200- 




200- 



Figure 15. MSSM. Left: Isolines of the total width T h3 (GeV). Right: isolines of BR(/i 3 -> hh). 
The light colored region is excluded at 95%C.L., the red region is excluded by CMS direct searches. 
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Figure 16. MSSM. Left: isolines of BR(/i 3 -> ri). Right: isolines of BR(/i 3 -> 66). The light colored 
region is excluded at 95%C.L., the red region is excluded by CMS direct searches. 



to be significant, especially for growing tan/3. In Figures 13 16 we draw the region excluded 



by such search, as inferred from 34 . 



6 The NMSSM at A > 1 and gauge coupling unification 

As said in the Introduction we are particularly interested in the NMSSM at A close to one 
and moderate tan /3 to limit the fine tuning. At least in the H decoupled case we have seen in 
Section [3] that this is consistent with current data. On the other hand a well known objection 
to the NMSSM at A > 1 is its compatibility with gauge coupling unification. Requiring A to 



stay semi-perturbative up to the GUT scale bounds A at the weak scale at about 0.7 35 . This 



value is in fact influenced by the presence of vector-like matter in full SU(5) multiplets that 
slows down the running of A by increasing the gauge couplings at high energies. However, 
even adding three vector-like five-plets at 1 TeV, in which case «g still remains perturbative, 



does not allow A at the weak scale to go above 0.8 36 37 



There are several ways 38-45 in which A could go to 1 -j- 1.5 without spoiling unification 



nor affecting the consequences at the weak scale of the NMSSM Lagrangian, as treated 
above. One further possibility makes also use of two vector-like five-plets as follows. For ease 
of exposition let us call them F u ^ + F Uj d, where F u is a 5 and Fj a 5, thus containing one 
S , f/(2)-doublet each, h u and hd, with the same quantum numbers of the standard H u ,Hd 
used so far. Correspondingly F U) d contain two doublets that we call h u ,d- Needless to say all 
these are superfields. Let us further assume that the superpotential is such that: 

• The five-plets interact with a singlet S and pick up S'C/(5)-invariant masses consistently 
with a Peccei-Quinn symmetry; 

• The standard doublets H u , Hd mix by mass terms with h u and hd, still maintaining 
the Peccei-Quinn symmetry, and do not interact directly with S. 
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The corresponding superpotential is 

/ = X s SF u F d + M U F U F U + M d F d F d + m u H u h u + m d H d h d + X t H u Qt, (6.1) 

where we have also made explicit the Yukawa coupling of the top to H u . Below these masses, 
all taken to be comparable, this / term leaves three massless supermultiplets: 

S, H u = c u H u + s u h u , H d = c d H d + s d h d , (6.2) 

which interact through the superpotential 

/ = \SH u H d + \ t H u Qt, A = \ s s u s d , X t = X t c u - (6.3) 

This superpotential, completed by Peccei-Quinn symmetry breaking terms at the Fermi 
scale, defines the effectiv e NM SSM as discussed so far. To cure the growth of A at increasing 
energies, the masses in (6.1) will have to be crossed while A is still semi-perturbative. For 
A = 1 -r 1.5 these masses are above 1000 TeV. 

At greater energies the running of the gauge couplings is affected, compared to the standard 
supersymmetric case, by the supermultiplet H u with a top Yukawa coupling increased by a 
factor l/c„ and by the degenerate complete S77(5)-multiplets F Utd + F Utd . To avoid a Landau 
pole before the GUT scale in the top Yukawa coupling c u has to be bigger than about l/\/2. 
As to the effect of F u>d + F u>d , they do not alter the relative one loop running of the gauge 
couplings but might give rise to an exceeding growth of all of them before Mqut due to the 
presence of the coupling As, which at some point will get strong. To avoid this a change 
of regime in the S*?7(5)-symmetric sector will have to intervene to keep under control the 
anomalous dimensions of the F u>d , F u>d superfields. 

7 Summary and conclusions 

Is the newly found Higgs boson at LHC alone or is it a member of a family of scalar particles? 
Now that we know that one Higgs boson exists this is a compelling question, both per se 
and in many motivated different contexts. Here we have focussed on the case implied by the 
supersymmetric extensions of the SM, where naturalness requires at least one extra doublet 
of Higgs particles at the Fermi scale. Specifically, a minimally fine tuned case appears to 
be the Higgs system of the NMSSM, which also includes an extra singlet state. It is in fact 
not unconceivable that part of such extended Higgs system be the lightest fragments of the 
entire new particle spectrum of the NMSSM, with the possible exception of the LSP. It is of 
particular interest to know which impact the measurements - current and foreseen - of the 
different signal strengths of the newly found resonance have on this problem and how they 
compare with the potential of the direct searches of the extra states. 

To have a simple characterization of the properties of the extended Higgs system we 
have focussed on relations between physical parameters. We have also not included any 
radiative effects of superpartners other than the top-stop loop corrections to the quartic 



Higgs couplings, (2.3), which mostly affect the properties of the lightest CP even state 
that we identify with /ilho While providing at least a useful reference case, we think that 
this is motivated by the consideration of superpartner masses at their "naturalness limit" . 
Furthermore we have considered two limiting cases at a time of the NMSSM, in which one CP 
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even scalar is decoupled: either the singlet or the doublet component, H = —spHd + c$H u . 



This allows to work in general with four effective parameters: A, tan/3, A t in (2.3) and the 
mass of the intermediate CP-even scalar, h 2 or h 3 in the two cases. The mass of hi = /ilhc se ts 
one relation among these parameters. A further relation exists in the singlet-decoupled case, 
of which the MSSM is a limit (A = 0). We have not been sticking to a particular NMSSM, 
which might imply specific constraints on the physical parameters that we consider, but we 
have assumed to live in the case of negligible CP violation in the Higgs sector. 

Within this framework, even though the signal strengths of /ilhc are close to those expected 
in the SM, they still allow for a new further state nearby, unlike in the case of the MSSM, 
where a CP- even scalar heavier than /ilhc and below about 300 GeV is unlikely [l~7] . This 
is true in both the limiting cases that we have considered, as visible in Figures [2] and |9l to 



be contrasted with Figure [T3j On the other hand, the same Figures show that the measured 
signal strengths of /zlhc do limit the possible values of A, at least for moderate values of tan /3, 
which is the region mostly motivated by naturalness: A ~ 1 is still largely allowed in the 
if -decoupled case, whereas it is borderline in the singlet decoupled situation. As commented 
upon in Section |6j we think that A > 1 can be compatible with gauge coupling unification. 
Most importantly from the point of view of the direct searches^ this feature reflects itself 
in the behavior of the new states, quite different in the two cases, especially in their decay 
properties. The state h 2 of the if -decoupled case has a large BR into a pair of h%, whenever 
allowed by phase space, with VV as subdominant decay (Figures [6]|5]). With the production 



cross sections shown in Figure K3^ its direct search at LHC8 or LHC14 may be challenging, 



although perhaps not impossible [47,48 . On the other hand the reduced value of A allowed 



in the singlet decoupled case makes the bb channel, and so the tt, most important, below 



the tt threshold (Figures 11 and 12). This makes the state h 3 relatively more similar to the 



CP-even H state of the MSSM (Figures 15 and 16), which is being actively searched. An 
important point in this comparison is that in the NMSSM there are two CP-odd states and 
that the lightest of the two may be quite different from the state A of the MSSM unless in 
a special S'-decoupled case. 

It will be interesting to follow the progression of the searches of the Higgs system of the 
NMSSM, directly or indirectly through the more precise measurements of the hi properties. 
We believe that the framework outlined here should allow to systematize these searches in a 
clear way. We also think that they should be pursued actively and independently from the 
searches of the superpartners. 
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